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ABSTRACT 

We present our analysis of UV attenuation by internal dust of a large sample (N=906 galaxies) 
of Lyman Break Galaxies (LBGs). Using spectral energy distributions (SEDs) from the PEGASE 
galaxy spectral evolution model we apply dust attenuation corrections to the G — 1Z colors using 
the Witt & Gordon (2000) models for radiative transfer in dusty galactic environments to arrive at 
the UV attenuation factors. We show that the dust in the LBGs exhibit SMC-like characteristics 
rather than MW-like, and that the dust geometry in these systems is most likely to be represented 
by a clumpy shell configuration. We show that the attenuation factor exhibits a pronounced 
dependence on the luminosity of the LBG, aigoo oc (X/Lq)", where 0.5 < a < 1.5. The exponent 
a depends on the initial parameters of the stellar population chosen to model the galaxies and 
the dust properties. We find that the luminosity weighted average attenuation factor is likely to 
be in the range from 5.7 — 18.5, which is consistent with the upper limits to the star formation 
rate at 2 < z < 4 set by the FIR background. This implies that the current UV/optical surveys 
do detect the bulk of the star formation during the epoch 2 < z < 4, but require substantial 
correction for internal dust attenuation. 

Subject headings: dust, extinction — galaxies: high-redshift — galaxies: ISM — ISM: evolution 



1. INTRODUCTION 

In the last few years, large samples of star form- 
ing galaxies have been observed at high redshift 
(z > 2) using the Lyman Break technique (Stei- 
del & Hamilton 1993). So far these surveys have 
been carried out mostly at optical wavelengths, 
which explores the rest-frame UV/FUV region of 
the spectral energy distribution (SED) of galaxies 
in the range 2 < z < 4. The UV luminosity in 
these star-forming galaxies relates directly to the 



number of short-lived, high-mass stars and this 
gives, in principle, a measure of the actual star 
formation rate (SFR). 

Attenuation due to internal dust presents one 
of the biggest obstacles to interpreting these ob- 
servations. Attenuation depends not only on the 
amount of dust, its composition and size distri- 
bution but also on its scattering properties and 
the geometry of the dust and its distribution rel- 
ative to the sources (Witt, Thronson, & Capuano 
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1992), none of which are well known for these 
high-z objects. The question of dust attenuation, 
therefore, becomes important when estimating the 
star formation rate history (SFRH) of the universe 
from observation of rest-frame UV fluxes of star- 
forming galaxies at high redshift. 

Many attempts have been made to estimate the 
SFRH in the recent years. One of the first was by 
Madau et al. (1996). This study, however, under- 
estimated the effects of dust. The fact that high 
redshift galaxies are intrinsically very dusty be- 
came apparent with the advent of the ISO obser- 
vations in the mid- and far-IR (Pugct ct al. 1996). 
Deep surveys in the FIR gave rise to galaxy counts, 
and COBE data on the FIR background resulted 
in upper limits for the SFR (Elbaz et al. 1999; 
Puget et al. 1999; Altieri et al. 1999; Aussel et al. 
1999). Even after correction for dust attenuation 
there still remained a discrepancy in the SFRs as 
determined by UV/optical surveys and the FIR 
background. There is a debate whether all the 
star-formation in the universe is in fact being de- 
tected by current UV/optical surveys. Poggianti 
& Wu (2000), Poggianti, Bressan, & Franceschini 
(2001) and Rigopoulou et al. (2000) report inde- 
pendent evidence on both local and high-z lumi- 
nous starbursts in which ~ 70% — 80% of the bolo- 
metric flux from young stars is completly obscured 
by dust and remains hidden in the UV/optical sur- 
veys (even after correction for dust). Adclberger 
& Steidel (2000) on the other hand suggest that 
the observed 850 fim galaxy counts and the back- 
ground can be explained by the LBG population 
by applying a proportionality correction to the op- 
tical flux, using the locally observed distribution 
of the mm-to-optical flux ratios. 

We believe these conflicts can be resolved by 
an improved knowledge of the dust, its properties 
and geometry in the high-z galaxies. We present 
our analysis for the attenuation correction for UV 
emission from a large sample (N=906) of Lyman 
Break Galaxies (LBGs), due to internal dust, at 
redshift 2 < z < 4. This is the largest sample 
of spectroscopically identified high-z galaxies to 
date. We rely on radiative transfer calculations 
to estimate the dust attenuation in each individ- 
ual galaxy. In this approach no assumptions are 
made to relate different SFR indicators, nor do we 
apply relations which have only been established 
in the local universe, to high-z objects. The only 



assumptions made are that these systems have a 
shell geometry and that the dust is clumpy, for 
both of which there are strong indications (Gor- 
don, Calzetti & Witt 1997). We explore the de- 
pendence of dust opacity on the intrinsic UV lu- 
minosity of the LBGs. We further investigate the 
possible evolution of the dust opacity in LBGs 
with respect to redshift and also the origin and 
nature of dust in the early universe. 

The layout of the paper is as follows. In section 
2 we describe the data and their limitations. Sec- 
tion 3 discusses the SEDs and the initial parame- 
ters used to model the intrinsic colors of the LBG 
population. In section 4 we describe the dust mod- 
els and the attenuation functions used. Section 5 
describes the analysis and the results obtained. In 
section 6 we present a detailed discussion compar- 
ing and contrasting our results with previous work 
as well as the implications of our results. Finally, 
section 7 summarizes our conclusions. 

2. DATA 

Our data consisted of 1Z- magnitudes, colors 
(U n — G,G — 1Z) and spectroscopic redshifts of 
~ 935 LBGs (Charles Steidel 2001, private com- 
munication) . The UV-dropout technique was used 
to identify these galaxies at z ~ 3 and spec- 
troscopic redshifts were obtained using the Low- 
resolution Imaging Spectrograph on the Keck tele- 
scopes. A detailed description of the filters and the 
observations can be found in Steidel & Hamilton 
(1993) and Steidel, Pettini, & Hamilton (1995). 
For rest-frame spectra of some of these objects see 
Pettini et al. (2001). Since the U-dropout was a 
selection criterion for the LBGs, the U n — G color, 
at best, was only a 1 sigma detection. Therefore 
we only used the G — 1Z color for our analysis. We 
also excluded the objects spectroscopically identi- 
fied as AGNs and QSOs, reducing the sample size 
to 906 galaxies. 

It was unfortunate that we were able to salvage 
only one color for such a large sample. Determin- 
ing the stellar content, dust content, dust prop- 
erties and the dust-star geometry, all from just 
one color is dangerous. Therefore, as a check for 
our model parameters, we applied our models to a 
smaller LBG sample for which longer colors were 
also available. We used the V and H magnitudes 
from Papovich, Dickinson, & Ferguson (2001) (31 
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galaxies) and Sawicki & Yee (1998) (17 galaxies). 
The H denotes the AB magnitudes from the NIC- 
MOS F160W filter and IRIM on the Kitt Peak 4.0 
m. V denotes the AB magnitudes from the HST 
F606W filter. The HST images were smoothed to 
match the much poorer PSF of the ground based 
IR observations. See Papovich, Dickinson, & Fer- 
guson (2001) and Sawicki & Yee (1998) for com- 
plete descriptions of the data. 

The G—1Z sample is subject to certain selection 
effects: only galaxies with G — 7Z < 1.2 were spec- 
troscopically studied and included in this sample, 
shown in Figure 1. This selection effect has sig- 
nificant consequences for the reddening estimates, 
particularly at large redshifts. No effects are in- 
troduced on the blue-limit of the observed colors, 
and the progressive reddening of the entire sample 
at z > 3 can be attributed to absorption by the 
intergalactic medium. (Madau 1995). 




2.8 3 3.2 3.4 
Redshift z 

Fig. 1. — G—1Z color plotted against spectroscopic 
redshift z for 906 LBGs. The upper cut-off at G — 
1Z ~ 1.2 is a selection limit. Progressive reddening 
at z > 3 is because of IGM absorption (Madau 
1995). 



3. SPECTRAL ENERGY DISTRIBU- 
TIONS 

The G — 1Z color of sources at z ~ 3 pro- 
vides a measure of the FUV spectral slope in their 
rest frame. The large range of G — 1Z colors in 
the LBGs at a given redshift suggests that the 
sources are reddened by different amounts and 
that this reddening is mainly due to dust intrinsic 
to the sources. It is therefore necessary to obtain 



an unreddencd spectral energy distribution (SED) 
from dust-free models of young star-forming pop- 
ulations for use as a template. We used the 
PEGASE2 galaxy spectral evolution model (Fioc 
& Rocca-Volmerange 1999) to model the LBG 
SEDs. The prominent feature of this model is 
that it includes stellar evolutionary tracks with 
non-solar metallicities and contribution from neb- 
ular emission. Zackrisson et al. (2001) provide 
a comparison of different galaxy spectral evolu- 
tion models used in recent studies. They show 
that the broad-band UV colors predicted by their 
model and the PEGASE2 codes are not signifi- 
cantly different for a large range of metallicities. 
There are differences in the V— K color, but this 
color is of little consequence in star-forming re- 
gions. After choosing the initial parameters to de- 
scribe an LBG, the resulting SED was redshifted 
by the measured redshift and convolved with the G 
and 72. transmission curves (Charles Steidel 2001, 
private communication; also see Steidel & Hamil- 
ton (1993)) to get the G — 1Z colors, using codes 
that are a part of the PEGASE2 package. 

The synthetic spectra are computed at specified 
ages after the onset of the star formation episode. 
The properties of the stellar populations are com- 
puted based on the shape of the initial mass func- 
tion (IMF) and its upper and lower limits. Other 
parameters used to compute the SED are the ini- 
tial metallicity of the ISM and the star formation 
history. 

Working with only one UV color, we face the 
dilemma produced by the age-reddening degener- 
acy. For a given source, we are unable to distin- 
guish between the cases of an older (and there- 
fore redder), dust-free population and that of a 
younger, dust-reddened population. There are 
various factors that can redden the colors of a 
galaxy; age of the stellar population, increased 
metallicity, higher amounts of dust or any com- 
bination of these factors. In principle, this de- 
generacy cannot be broken working with one UV 
color alone. In our analysis we break this degen- 
eracy arbitrarily by analyzing the sample under 
two sets of assumptions regarding the time depen- 
dence of dust formation in these high-z galaxies, 
as discussed in the next two subsections. In each 
of the two cases we assume a fixed age of the stel- 
lar population for all the LBGs and assume con- 
ditions that leads to an upper-limit and a lower 



3 



limit on the dust attenuation. We find that only 
a small fraction of the color range of the LBGs 
can be attributed to age alone; most of it must be 
due to reddening by internal dust. As shown in 
Figure 2, at a given rcdshift (z <~ 2), the G — 1Z 
color gets redder by 0.26 as the galaxy ages from 
yr to 1 Gyr and by 0.15 as the galaxy ages from 
50 Myr to 1 Gyr. The color range for the ob- 
served LBGs at the same redshift is 1.6. Analysis 
under two different sets of assumptions based on 
the timescale of dust formation helps in resolving 
this age-reddening degeneracy as far as it exists. 
Furthermore we verify our model parameter as- 
sumptions by applying them to a subset of the 
LBGs for which longer colors (V — H) were avail- 
able. At z > 3, the colors of a galaxy arc fur- 
ther reddened by IGM absorption. Both optically 
thin Lyman— a forest clouds and optically thick 
Lyman-limit systems contribute to the reddening 
of background sources (Madau 1995). We have 
corrected our assumed SEDs for this effect using 
the values provided in Madau (1995). 
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Fig. 2. — G — 1Z color predicted by increasingly 
older SEDs as a function of redshift z. The range 
on observed colors of the LBGs cannot be ex- 
plained by age alone. 

3.1. Case I 

In this case we assumed that we are observ- 
ing the LBGs within a short time into their first 
episode of star formation. The dust free model 
was chosen to represent the very youngest, bluest 
LBG. This implies that only the very bluest of 
the LBGs are dust free and that increased red- 



dening is a consequence of rapidly forming dust. 
We further assumed that dust formation follows 
the formation of massive stars on a nuclear time 
scale. In order to get an absolute lower envelope 
to the data at higher redshifts (z > 2.5) we had 
to assume lower metallicities (and correspondingly 
younger ages). The initial metallicity of the ISM 
was varied from 0.1 Z Q to 0.01 Z Q with z = 2 to 
z = 4 respectively. In order to mimic the behavior 
of the data at the blue envelope we had to in- 
crease the higher mass limit for the IMF at larger 
redshifts. We used a Salpeter initial mass function 
with a lower mass limit of 0.1 M@. Studies of the 
early metal pollution of the intergalactic medium 
show that the IMFs are increasingly biased to- 
ward higher masses at large redshifts (see Abia 
et al. (2001) and references therein). Also Lloyd- 
Ronning, Fryer, & Ramirez-Ruiz (2002) show that 
if massive stars are the progenitors of long lived 
gamma ray bursts, the IMFs at large redshifts 
need to be flatter with much higher upper mass 
limits. The upper mass limit for our models was 
varied from 120 M Q to 200 M Q with z = 2 to 
z = 4 respectively. We divided the sample into 
four redshift bins (2 < z < 2.5, 2.5 < z < 
3, 3 < z < 3.5, and 3.5 < z < 4), to simulate 
the changing behavior of the data with redshift as 
well as to study the possible evolution of the dust 
properties. Table 1 shows the parameters which 
were required to match the bluest observed col- 
ors with an unreddened SED for case I. A con- 
stant Star Formation Rate (SFR) of 100 M Q y^ 1 
was used at all redshifts. The colors were exam- 
ined at different ages after the star formation had 
started as noted in Table 1 under the column la- 
beled "Age". Values for the UV-spectral slope /3 
were calculated as prescribed by Calzetti (2001) 
for the spectral range 0.126 [im < A < 0.26 fim. 
A scenario favoring an IMF biased toward mas- 
sive stars, low stellar metallicity, and substantial 
dust formation on a 10 yr time scale is supported 
by the equivalent width distribution and number 
counts of Lyman— a emission galaxies in the early 
universe by Malhotra & Rhoads (2002). The in- 
trinsic colors of all the galaxies are assumed to 
correspond to this young model. 

In this case the model parameters were cho- 
sen so as to provide a blue envelope to the sam- 
ple, representing the youngest unreddened galax- 
ies. In doing this, we have assumed a short time- 
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scale for dust formation and assumed that only 
the youngest systems are dust free. In such young 
star forming systems Type II supernovae (SNe II) 
are the dominant sources of dust production. (For 
a detailed study on the observed properties of SN 
dust see Wooden (1997)). Modeling the dust con- 
tent of a nearby low mctallicity galaxy incorpo- 
rating dust from SNe II Hirashita, Hunt, & Fer- 
rara (2002) show that the dust mass reaches its 
maximum in ~ 2 x 10 7 yr. Edmunds & Eales 
(1998) model dust masses in galaxies and show 
that the dust content in galaxies with high star 
formation rates advances within a factor of 2 of 
the maximum values at very young ages. For 
m+/mtot = 0.2, if one considers m to t = 10 11 Mq, 
and SFR = 100 M^yr' 1 , dust masses grow within 
a factor of 2 of the maximum in 2 x 10 7 yr. This 
assumption of the short time scale for dust forma- 
tion results in larger dust attenuation factors for 
these galaxies. 

3.2. Case II 

In this case we assumed longer timcscales for 
dust formation, consistent with the assumption 
that mass loss from evolved lower-mass stars is 
the principal contributor to dust formation. This 
process follows the slower (r ~ 10 9 yr) evolution- 
ary time scale of the stars; it is more significant in 
the present universe and less so in the high-z uni- 
verse (Edmunds 2001). The model parameters for 
this case were chosen, so as to arrive at a conser- 
vative lower limit to the amount of dust present in 
LBGs responsible for partial attenuation of their 
UV flux. Shapley et al. (2001) analyzed 81 LBGs 
for which they also had K s and J band fluxes, and 
these were found to have ages between 10 Myr and 
1 Gyr with a median value of 320 Myr. We as- 
sume that galaxies with populations younger than 
320 Myr are dust free and that dust appears only 
through subsequent evolution. The 320 Myr old 
SED is therefore used as a reference for all LBGs 
with colors redder than those implied by this ref- 
erence. As we were not attempting to model the 
bluest LBGs, dividing the sample into different 
rcdshift bins would not constrain the model pa- 
rameters in any way. Therefore, in this case we do 
not divide the sample into different redshift bins. 
A constant SFR of 100 M Q yr^ 1 and a Salpeter 
IMF with lower limit of 0.1 Mq and an upper limit 
of 125 Mq were chosen, which are consistent with 



those found by Shapley et al. (2001). The initial 
metallicity of the ISM was set to 0.2 Zq and the 
colors were sampled at 320 Myr. This larger age 
and higher metallicity vis-a-vis case I result in red- 
der colors for the dust-free (unreddened) galaxies, 
and all the galaxies with colors bluer than this 
were assumed to be dust-free. Table 1 shows the 
parameters used for the analysis of this case. 

For this case, we assume that the galaxies 
320 Myr old and younger are still dust free. Dust 
forms in these galaxies at later times. This long 
time-scale for the dust formation implies that the 
dust formation in these galaxies was primarily a 
result of the evolved intermediate mass stars and 
the dust production time-scales are linked to the 
evolution time-scales of intermediate mass stars. 
As explained earlier, we adopt an age of 320 Myr 
to account for the fact that Shapley et al. (2001) 
found that this was the median age of LBGs us- 
ing a multi-wavelength study, and this approach 
results in a lower estimate to the amount of dust 
in these systems. 

These two cases are not just two possible sce- 
narios, but from the viewpoint of dust forma- 
tion scenarios are extreme models, and their range 
of results will cover any reasonable intermediate 
model. 

4. DUST MODELS AND ATTENUA- 
TION FUNCTIONS 

To relate the SEDs generated with the assump- 
tions outlined in Sec. 3 to the observed colors, we 
need to know the wavelength dependence of the 
dust attenuation function in these galaxies. It is 
important to distinguish between the terms ex- 
tinction and attenuation. The term "extinction 
law " was originally defined for stars. It quanti- 
fies the wavelength dependence of dust absorption 
and scattering out of the line of sight toward point 
sources. In distant galaxies however the situation 
is quite different. There are two main reasons: (i) 
The complex geometry of the light sources and the 
dust distribution and (ii) the flux scattered off the 
dust grains is returned into the line of sight in com- 
plex ways. Thus the "attenuation" depends criti- 
cally on the geometry of the dust and the sources 
and the optical depth and scattering properties of 
the dust. 

In recent years the wavelength dependence of 
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Table 1 

Parameters assumed for SEDs under Case I and Case II 



Rcdshift Range Initial Metallicity Age Upper mass limit (3 



Case I 
2.0 < z < 2.5 
2.5 < z < 3.0 
3.0 < z < 3.5 
3.5 < z < 4.0 

Case II 
2.0 < z < 4.0 



0.1 Z 
0.1 Z Q 
0.01 Z Q 
0.01 Z Q 



0.2 Zr. 



- 50 Myr 

- 50 Myr 
~ yr 
~ yr 

320 Myr 



120 Mq 
135 M 

150 Mq 

200 M 



125 Mc 



© 



-3.14 ±0.04 
-3.5 ±0.02 
-3.76 ±0.02 
-3.82 ±0.02 



-3.04 ±0.04 



dust attenuation in star forming galaxies has been 
estimated by the so called "Calzetti law", which 
is an empirical determination of the observed at- 
tenuation function derived by averaging over a 
small sample of UV-bright local starburst galax- 
ies (Calzetti, Kinney & Storchi-Bergmann 1994; 
Calzetti et al. 1995, 2000). This empirical law 
leads to a correlation between the UV spectral 
slope (3 and the FUV attenuation. The FUV flux 
has then been corrected for attenuation, assum- 
ing this correlation holds for all redshifts. This 
law however has limitations; while applicable to 
the particular sample there is no evidence that 
this law applies over all redshifts or that it applies 
to systems with substantially more dust than is 
present in the local sample of UV-bright galaxies. 
Bell (2002) demonstrates that normal star-forming 
galaxies deviate substantially from the starburst 
(3— Apuv correlation. We have therefore chosen to 
address the question of dust attenuation with the 
radiative transfer models of Witt & Gordon (2000) 
(henceforth WG2000). These models incorporate 
multiple scattering radiative transfer calculations 
for different galactic environments, filled with ei- 
ther homogeneous dust or a two-phase clumpy 
dust distribution. Figure 3 shows the attenuation 
functions using the clumpy shell models for various 
optical depths. In these models, ry should be in- 
terpreted as a measure of the dust mass (ry oc dust 
mass) where the actual amount of dust would scale 
with the size of the system. As the amount of dust, 
measured by the visual optical depth for a corre- 
sponding homogeneous model increases, the atten- 



uation law becomes flatter in the far-UV. Thus, as 
Figure 3 shows, there is not a single attenuation 
law, e.g. the Calzetti "law", but a series of at- 
tenuation laws with a functional form dependent 
upon the amount of dust present. WG2000 have 
also shown that the Calzetti law corresponds to 
the attenuation function for a clumpy shell galaxy 
with SMC dust of intermediate dust optical depth 
of ry = 1.5 (see Figure 3). One cannot expect all 
the LBGs to have the same amount of dust given 
their wide range of observed colors. 

The clumpy models of WG2000 assume a spa- 
tial filling factor of 0.15 for denser clumps, em- 
bedded into a low density inter-clump medium, 
with a density contrast of 100 : 1. As shown 
in an earlier study (Witt & Gordon 1996), this 
set of conditions leads to a power-law spectrum 
of cloud masses that closely resembles that found 
in the galactic ISM. The dependence of radiative 
transfer on other choices of dumpiness parameters 
was fully explored by Witt & Gordon (1996). In 
the absence of any other information, we assume 
that the structural details of the ISM in LBGs 
resemble that of the MW galaxy. The geometry 
of the WG2000 shell models assume an inner vol- 
ume (r < 0.3R) occupied by stars, surrounded by 
a shell (0.3R < r <R) of dust distributed either 
homogeneously or in a two-phase clumpy distri- 
bution. While relative sizes of star-filled volume 
and the dusty shell are to some degree arbitrary, 
they are modeled after he appearance of giant star- 
forming complexes, e.g. NGC 604 in M33 (Hunter 
et al. 1996). Important aspects are that the dust 
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Fig. 3. — Attenuation functions from clumpy shell 
models at different optical depths ry compared to 
the Calzetti law. Dotted line is the attenuation de- 
scribed by E(X - V)/E(B - V) for MW-like dust 
at Ty — 1.5. Long-dashed line is the attenuation 
curve for SMC-like dust at ry = 0.5, short-dashed 
for SMC-like dust at ry = 1.5, and dash-dotted for 
SMC-like dust at ry = 4.5. Solid line represents 
the Calzetti law where the solid squares are photo- 
metric points. The Calzetti law can be reproduced 
by SMC-Clumpy shell model with ry = 1.5. 



is situated between the sources and the observer, 
thus maximizing the attenuation, yet close enough 
to the sources so that scattered light becomes an 
integral part of the flux received from the complex. 

In this study we have used the clumpy shell 
models to estimate the dust attenuation for the 
LBGs (See sec. 5.1 for a more detailed analysis 
of the dust geometry). Figure 4 shows the rela- 
tionship between the attenuation at 1600 A and 
the optical depth for SMC-clumpy shell and SMC 
- homogeneous shell models. In the homoge- 
neous models the dust has uniform spatial density. 
This effectively represents a screen geometry, and 
the attenuation rises linearly with optical depth. 
Whereas in the two-phase clumpy models, at low 
values of ry both the dust clumps and the inter- 
clump medium are optically thin and the attenua- 
tion increases proportional to the amount of dust, 
similar to the homogeneous case. At higher val- 
ues of ry the clumps are starting to become opti- 
cally thick but the inter-clump medium is still op- 
tically thin and the attenuation varies less rapidly 
with the amount of dust. At intermediate optical 
depths the attenuation is influenced both by the 
dust clumps and the inter-clump medium. 

5. ANALYSIS 

5.1. The Nature and Geometry of the 
Dust in LBGs 

There are significant variations in the dust 
properties in different galactic environments of co- 
eval galaxies within the local group itself. In the 
LMC, the UV extinction curves show a distinctly 
different behavior between the 30 Dor region (a 
mini-starburst) (Walborn 1991) and the rest of the 
LMC (Clayton & Martin 1985; Fitzpatrick 1985, 
1986; Misselt, Clayton, & Gordon 1999). The 2175 
A bump is weaker and the far-UV rise is stronger 
near the 30 Dor region than in the rest of the LMC. 
In the bar-region of the SMC, the average extinc- 
tion curve is characterized by a roughly linear rise 
(vs. A" 1 ) increasing toward shorter wavelengths 
without a 2175 A bump (Prevot et al. 1984; 
Thompson et al. 1988; Gordon & Clayton 1998). 
Yet, there is one sight line that has an extinction 
curve with a significant 2175 A bump (Lequeux 
et al. 1982; Gordon & Clayton 1998). In M31, 
the extinction curve is consistent with that of the 
average Galactic extinction within the associated 
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uncertainties, although the 2175 A bump may be 
weak (Bianchi et al. 1996). Using IUE data for 
sight lines through low-density regions of the MW 
disk and halo Clayton, Gordon, & Wolff (2000) 
show that many of the sight lines have extinction 
curves with weak bumps and very steep far-UV 
extinction reminiscent of the Magellanic Clouds. 
Modeling the SEDs of 30 starburst galaxies Gor- 
don, Calzetti & Witt (1997) found that the dust 
in these systems has extinction properties similar 
to the SMC lacking the 2175 A feature present 
in the galactic extinction curve. However, Motta 
et al. (2002) report on the detection of a signif- 
icant 2175 A feature in the extinction curve of 
a normal-type galaxy at z = 0.83. Massarotti et 
al. (2001a) report that some galaxies in the HDF 
spectroscopic sample show the 2175 A feature. 

The limited evidence available so far suggests 
that the absence of the 2175 A feature is associ- 
ated with high SFR, but this is far from conclusive. 
In the presence of this ongoing debate, the nature 
of the dust in the LBGs needs to be addressed. We 
use the fact that the presence of the 2175 A fea- 
ture has a drastically different effect on the G — TZ 
color in the 2 < z < 2.6 range compared to dust 
without the 2175 A feature to test the LBG data 
for the presence of the 2175 A band. We red- 
dened the assumed SED (case I) with the WG2000 
clumpy attenuation functions with increasing op- 
tical depths of t v = 0.25, 0.5, 1.0, 5.0, 7.0, & 10.0, 
with both MW- and SMC-like dust. These SEDs 
were then redshifted and the G—TZ color was calcu- 
lated at the different rcdshifts. For z > 3 the col- 
ors were corrected for IGM absorption using values 
in Madau (1995). Figure 5 is a plot of the G - TZ 
color vs. the redshift, where the different lines are 
the colors obtained for the intrinsic SED reddened 
with different amounts of MW-like dust. We find 
that MW-like dust cannot explain the range of ob- 
served colors in the 2 < z < 3 range. As the 
2175 A feature passes through the G and TZ fil- 
ter set at different z, the predicted G — TZ color 
first gets bluer (till z < 2.45) and then gets pro- 
gressively redder for increasing amounts of dust 
(see Figure 5). On the other hand, the curves 
modeled by SMC-like dust (Fig. 6) for increas- 
ing values of attenuation span the observed range 
of colors quite well. Using color-color plots Gor- 
don, Smith, & Clayton (1999) also reached the 
conclusion that the dust in starburst galaxies in 



the North and South Hubble Deep Field is similar 
to SMC-like dust. Figure 6(a) shows the model- 
predicted curves for SMC clumpy dust under the 
case I assumptions for the intrinsic SED, and Fig- 
ure 6(b) for the case II assumptions. We conclude 
that the dust in this sample of LBGs is better rep- 
resented by SMC-like dust, which lacks the 2175 
A feature. 

Even though multi-color information is avail- 
able for some LBGs, such samples are much 
smaller and therefore less complete, being ad- 
versely affected by observational limitations. We 
have chosen to work with the largest sample of 
LBGs available to date. Unfortunately this im- 
plies that we infer dust content, dust properties, 
stellar content and dust-star geometry all with 
a single color. As a check to our model SED 
assumptions, we analyzed a smaller sample of 
LBGs for which longer wavelength colors were 
available. Figure 7 is the V — H color vs. red- 
shift for LBGs (data from Papovich, Dickinson, & 
Ferguson (2001) and Sawicki & Yee (1998)). The 
curves on the figure are model-predicted curves 
for SMC clumpy dust, with the dust-free intrin- 
sic SED produced under the same assumptions 
as our case I analysis. The colors predicted by a 
dust-free SED under case II assumptions is also 
shown for reference. Comparison of Figure 6(a), 
6(b) with Figure 7 clearly shows that the smaller 
V — H sample spans the same range in attenua- 
tions and that our model parameters are robust. 
This also implies that these two samples of LBGs 
are subsamples of galaxies in the same evolution- 
ary stage. 

We need a fundamental understanding of the 
environment that produces the large range of red- 
dening that is observed. The dust geometry is a 
important factor in the amount of reddening pro- 
duced. For a local sample of starburst galaxies 
Gordon, Calzetti & Witt (1997) have shown that 
the dust geometry is best represented by an in- 
ner dust-free sphere of stars surrounded by a star- 
free shell of clumpy dust such as represented by 
the WG2000 shell models. Inoue (2002) has found 
the dust in HII regions to be best described by 
a shell geometry. Also Buat et al. (2002) find 
that the clumpy shell geometry as described in 
WG2000 fits very well with their analysis of the 
UV SEDs of starburst galaxies. Different dust ge- 
ometries are analyzed in WG2000. Not all geome- 
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f SMC - Clumpy Shell '■ 

/ SMC - Homogeneous Shell ■ 




Fig. 4. — Attenuation at 1600 A , A1600 ( ma g) 
against dust optical depth, Ty for different ge- 
ometries. Solid line represents SMC-like dust in 
a Clumpy-shell geometry and dashed line repre- 
sents SMC-like dust homogeneously distributed in 
a shell geometry. 




Redshift z 

Fig. 5. — G — 1Z color vs. redshift where the 
different lines are the colors obtained for the 
intrinsic SED reddened with different amounts 
of MW-like dust. The attenuation (in mag- 
nitudes) produced at 1600 A , ^4i6oo is indi- 
cated for each of the lines. The values of Ty 
in the model that produce these attenuations are 
0, 0.25, 0.5, 1.0, 5.0, and 10.0. Even large amounts 
of dust cannot explain the red colors in the 2 < 
z < 3 region. 




Fig. 6. — G — 1Z color vs. redshift with the col- 
ors obtained for the SED reddened with differ- 
ent amounts of SMC-like dust, (a) under Case I 
assumptions and (b) under Case II assumptions. 
The attenuation produced at 1600 A is indi- 
cated for each of the lines. The values of Ty in 
the models that produce these attenuations are 
0,0.25, 0.5, 1.0, 5.0, 7.0, and 10.0. 
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tries produce the same amounts of reddening for 
given amounts of dust; to produce equal amount 
of reddening other geometries would require much 
larger dust masses. In their Figure 11 WG2000 
show that only the shell models reproduce as large 
a range of A/3 as observed. One should note that 
the observed G — 1Z color and (3 are not identi- 
cal. The G — 7Z color is a rest-frame UV-color, 
whose rest-frame UV base wavelengths vary with 
redshift (G = 1610 A at z = 2 to 966 A at z = 4, 
K = 2310 A at z = 2 to 1386 A at z = 4), while 
the wavelengths for (3 are fixed (1200 A-2600 A). 

It is not surprising that a large fraction of the 
galaxies at z > 2 are best described as resem- 
bling centrally concentrated "blobs", which arc 
probably the progenitors of our present-day galac- 
tic bulges and elliptical galaxies (van den Bergh 
2002). The choice of a spherical star-forming 
region, relatively free of dust surrounded by a 
dust shell seems to be an appropriate geometry 
to describe these objects. A justification for the 
clumpy-shell geometry of the LBGs is due to their 
large SFRs, up to ~ 1000 M yr _1 in some cases. 
Such massive star formation activity will blow 
away the gas and dust into the surrounding region. 
Stars are expected to spend the very beginning 
of their evolution deeply embedded in dusty envi- 
ronments, later dispersing the molecular clouds in 
which they were born (Calzetti, Kinney & Storchi- 
Bergmann 1994, and references therein). In our 
analysis too, we find that the shell geometry best 
explains the observed colors for reasonable values 
of dust optical depths. 

5.2. The Relation Between the Dust Opac- 
ity and the Intrinsic UV Luminosity 
of LBGs 

To estimate the intrinsic UV luminosity of the 
LBGs we need to adopt a cosmology: throughout 
this paper we use H = 65 km s" 1 Mpc^ 1 , Qm = 
0.3, f^A = 0.7, unless otherwise specified. The red- 
dened SEDs were redshiftcd and the colors were 
calculated for different amounts of dust, result- 
ing in a range of attenuations at 1600 A (Aiqqq), 
at different redshifts. The attenuation for each 
galaxy was then computed. This was done as fol- 
lows: The color curves corresponding to selected 
values of the attenuation at 1600 A , Ai e00l are 
as shown in Figure 6(a) and 6(b). Dividing the 
sample into finite-sized ranges of color surround- 



ing each curve partially accounts for color differ- 
ences which may be attributed to a range in the 
intrinsic ages of the SED. The attenuation for each 
galaxy in the rest-frame wavelength for the 1Z filter 
is then noted. The attenuations and z are used to 
arrive at the absolute magnitudes, luminosity dis- 
tances and the intrinsic UV luminosities. We find 
that the bulk of the LBGs are extremely luminous, 
L ~ 10 11 - 1O 12 L , corresponding to the LIRG 
and ULIRG population in the present universe. 

Star forming galaxies, both in the local universe 
and in the early universe have been shown to ex- 
hibit a correlation between dust attenuation and 
their intrinsic luminosity. Adelberger & Steidel 
(2000), in their Figure 11, relate the dust lumi- 
nosity to the total luminosity of starbursts in the 
local universe, at z <~ 1 and at z ~ 3. Calzetti 
et al. (1995), Sullivan et al. (2001), and Wang 
& Heckman (1996) report similar correlations be- 
tween different star formation indicators to the 
FIR luminosity (which is directly related to the 
amount of dust). Hopkins et al. (2001) analyze 
the effect of dust- reddening dependent on SFR ap- 
plied to different SFR indicators such as UV lumi- 
nosity, H a flux, FIR Luminosity and radio lumi- 
nosity at 1.4 GHz. They show that discrepancies 
between these methods can be explained by such 
an effect. All these trends imply that the more 
luminous galaxies are dustier. 

In this study also we find a correlation between 
the intrinsic UV luminosity and the dust atten- 
uation in these systems. The LBGs have almost 
constant apparent magnitude (~ 24.5) but a wide 
range of colors. This seems to imply a correlation 
between the intrinsic luminosity of these galaxies 
and their reddening. Given our large sample and 
systematic analysis we can arrive at a quantitative 
result, but we should point out that this relation 
is not unique; it is dependent upon the dust ge- 
ometry, dust properties and the assumed intrinsic 
SED. We find, 

~- f ^isoo A \ / 1 \ 

ai600 oc I — — — I (1) 

where ai6oo is the linear attenuation factor re- 
lated to the attenuation Ai 6 oo (in mag) as ai 6 oo = 
10°- 4j4l60 °. The value of a is slightly different for 
case I and case II. For case I, the mean rest-frame 
UV luminosity of all the galaxies in each attenua- 
tion interval was calculated for each redshift bin. 
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A similar procedure was followed for case II, ex- 
cept that there was no redshift binning. Figure 8 is 
a plot of Li6oo/£© vs. the attenuation at 1600 A, 
for case I. We find that the attenuation factor and 
the intrinsic UV luminosity are correlated, with 
a = 0.90 ±0.55 

This does not imply that the dust masses in- 
crease nearly linearly with UV luminosity. In 
a clumpy shell geometry the actual dust masses 
(measured by ry) increases (from ry = 1 to 
Ty = 10) by a factor of 10 while the attenuation 
increases (from 1.8 to 4.6 mag) by a factor of 13.2 
(2.8 mag). As Figure 4 shows, the Ai &m — T v 
relation is model dependent, and the actual dust 
masses can be inferred from the observed attenu- 
ation for a particular geometry. In homogeneous 
shell geometries the attenuations would increase 
much more rapidly for a small increase in the dust 
masses. 

As is shown in Figure 8, there is no significant 
evolution of the attenuation-luminosity relation- 
ship with redshift; the relations in the four red- 
shift ranges are essentially the same within the 
scatter of the data. This implies that the dust 
is mainly a direct by-product of the currently on- 
going star formation. As the star-forming region 
ages, new stars are formed with increasing initial 
metallicities, reflecting the increasing metallicity 
of the ISM. However, the dust is mainly produced 
by the massive stars and SNe II and the rate of 
dust production or the amounts of dust produced 
does not evolve significantly (see sec. 3.1 and sec. 
3.2 for a detailed discussion of the origin of the 
dust in LBGs). It should however be noted that 
because of the selection effects at the red end of 
the G — 1Z color range, both the 3.0 < z < 3.5 and 
3.5 < z < 4 bins are incomplete and the relation- 
ship may evolve beyond z ~ 3.0. 

Figure 9 is the plot of Liq 00 /L q vs. the atten- 
uation at 1600 A, for case II. Here one SED was 
assumed to represent the intrinsic SED of the en- 
tire population, with a longer timcscale for dust 
evolution. As can be seen from Figure 6(b), many 
more galaxies are assumed to be dust-free than in 
case I. This changes the dependence of the intrin- 
sic luminosity on the dust opacity. For this case 
we find a = 0.97 ±0.49. 
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Fig. 7. — V — H color vs. redshift where the differ- 
ent lines are the colors obtained for an SED red- 
dened with different amounts of SMC-like dust. 
The data points denoted by 'plus' symbols are 
from Papovich, Dickinson, & Ferguson (2001) and 
those denoted by solid circles are from Sawicki & 
Yee (1998). The line styles and symbols used are 
the same as those used in Figure 5 and Figure 6. 
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Fig. 8. — Mean Luminosities (ii6oo/-^©) m differ- 
ent attenuation (Aiqqq) intervals as calculated un- 
der Case I. The line is the mean of the least-square 
fits for the luminosity-attenuation correlation, for 
the different redshift ranges, with a = 0.90 ±0.55. 
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5.3. Luminosity Weighted Attenuation 
Distribution Functions 



Given the correlation between dust opacity in 
the UV and the intrinsic UV luminosity of LBGs 
and the wide range of derived attenuation val- 
ues, the appropriate way to present the average 
attenuation factor in these galaxies is a luminos- 
ity weighted average attenuation factor. A nor- 
malized luminosity weighted attenuation distribu- 
tion function for these galaxies is instructive. Fig- 
ure 10 is the distribution function of attenuations 
for LBGs under case I. As no redshift dependence 
was found in our analysis, the weighted attenua- 
tions for each attenuation interval were added at 
all redshifts. As the Figure clearly illustrates there 
is a large range of attenuations that these galaxies 
suffer, with the bulk of the luminosity attenuated 
by a factor of 10 - 50. The assumption that all 
the observed galaxies are very young and are be- 
ing observed during their very first star-formation 
episode, leads to the interpretation that all the 
reddening observed is a result of dust attenuation. 
Thus, the average value of the luminosity weighted 
attenuation factor for the whole sample, 18.4 (3.17 
mag), is considered to be an upper limit for this 
sample. 

Figure 11 is the distribution function of attenu- 
ations under case II. This sample, also, exhibits a 
large range in attenuations, with the bulk of the lu- 
minosity however suffering lower attenuations and 
only a small population being heavily attenuated. 
The average value of the luminosity weighted at- 
tenuation factor under this case is 5.7 (1.88 mag). 
This value can be understood as a lower limit of 
the attenuation factor for the sample as a whole. A 
comparison of these two figures clearly illustrates 
the differences which are a result of the different 
initial assumptions. 

6. DISCUSSION 

6.1. The UV Attenuation Factor and its 
Implications for the Star Formation 
Rate History of the Universe 

Our results have significant implications for the 
star formation history of the universe. UV lumi- 
nosity in star-forming galaxies relates directly to 
the number of short-lived, high-mass O & B stars, 
and this gives, in principle, a measure of the actual 
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Fig. 9. — Mean Luminosities (L1600/L©) in dif- 
ferent attenuation (A1600) intervals as calculated 
under Case II. The line is the least-square fit to 
the data, with a = 0.97 ± 0.49. 
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Fig. 10. — Normalized luminosity weighted atten- 
uation distribution function under case I. The bulk 
of the galaxies are attenuated by 2.5-4.25 mag 
(factor of 10-50). 
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Fig. 11. — Normalized luminosity weighted atten- 
uation distribution function under case II. The 
bulk of the galaxies are attenuated by smaller 
amounts. 

SFR after appropriate correction for dust attenu- 
ation. Therefore, the estimates of SFR at high-z 
are closely related to the FUV attenuation factor. 

In recent years many estimates for the FUV at- 
tenuation factor at high-z have been made based 
upon different methods. Based on metallicity ar- 
guments Madau et al. (1996) estimated the atten- 
uation factor to be zero. Very soon however it 
was recognized that the UV luminosity is substan- 
tially suppressed by dust, and using the Calzetti 
law to estimate the attenuation factor from the ob- 
served UV-slope, many have found values ^4—15 
(Meurer et al. 1997; Meurer, Heckman, & Calzetti 
1999; Steidel et al. 1999; Massarotti, Iovino, & 
Buzzoni 2001b). All these authors exploit the lo- 
cally calibrated relationship between the ratio of 
FIR and UV fluxes and the UV spectral slope /3. 
As shown in section 4 this law can be reproduced 
by models for galaxies with intermediate optical 
depths, but there is no reason to believe that all 
the galaxies have similar UV-slopes or an unique 
value of the dust optical depth, given the wide 
range of colors they exhibit. Pettini et al. (1998) 
obtained an attenuation factor of 2.5 — 6.3 for a 
sample of 5 LBGs, by estimating SFR using 
luminosity and requiring that SFR^ =SFR' [/y , 
where SFR' are the SFRs corrected for attenua- 
tion. Estimates of the SFR derived from the EFj lu- 
minosity are corrected for dust attenuation based 
upon the change in the Balmer decrement. These 
lines are produced in the ionized envelopes of O 



Stars and not B stars. Furthermore the stellar 
continuum, which would be measured in SFR'^y, 
suffers only about half the reddening suffered by 
the ionized gas (Calzetti 2001), which makes this 
method unreliable. Nandra et al (2002) and Seib- 
ert, Heckman, & Meurer (2002) use the X-ray lu- 
minosity to derive UV attenuations, using the as- 
sumption that the locally derived correlations be- 
tween X-ray luminosity and bolometric luminosity 
holds for all rcdshifts. This is unlikely because the 
sources of X-ray luminosity evolve through time; 
in the local universe the X-ray luminosity is domi- 
nated by low mass X-ray binaries whereas at large 
redshifts supernova remnants of massive stars are 
likely to be the dominant X-ray sources. See Ta- 
ble 2 for different values for the UV attenuation 
obtained in the recent years. 

Most of these methods arrive at the UV atten- 
uation from observations at different wavelength 
regimes and then use the Calzetti attenuation 
curve to arrive at the UV attenuation values. Use 
of a single attenuation curve for all the galaxies 
at all redshifts is not valid. Moreover, estimates 
based on small samples probably do not represent 
the whole. Working with larger samples has given 
rise to a large range in these estimates even by 
these methods (Meurer et al. 1997; Massarotti, 
Iovino, & Buzzoni 2001b). Also, these authors 
quote simple average values of the attenuations 
whereas we feel that luminosity- weighted attenua- 
tion values are more valid as the star formation in 
these high- z systems occurs in galaxies with a wide 
range of luminosities and is apparently dominated 
by star- formation in extremely luminous systems. 

Quite encouragingly, our results of dust attenu- 
ation correction of the UV data is consistent with 
the constraints imposed by the cosmic infrared 
background (CIRB). Chary & Elbaz (2001) use a 
variety of evolutionary models to fit the local mid- 
IR luminosity function of galaxies with the CIRB. 
They deduce that the dust enshrouded SFR den- 
sity peaks at z — 0.8 and remains flat up to a red- 
shift of z ~ 4. In Figure 12 the three curves are 
the SFRs predicted by three FIR models (Chary 
& Elbaz 2001) compared to the observed SFR 1 
(Steidel et al. 1999) and the attenuation corrected 
values. The lower limit on our points is obtained 



1 Cosmology used to estimate SFR: Qm 
0.7, H = 50 km s" 1 Mpc" 1 



0.3, n A 
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Table 2 

Different values for the UV attenuation factors obtained in recent years 



Method 


KJ Cll 11 L7 il)lzj(_< 


TTV ^ttoTinstion factor 


A Ivjlvl V. 11 V.. A. 


Metallicity 


60 





1 


Calzetti law 


60 


3-15 


2 


Ha luminosity 


5 


2.5-6.3 


3 


Calzetti law 


25 


5.4 


4 


UV slope 


48 


4.8 


5 


Multi-wavelength analysis 




8 


6 


III emission 


5 


3 


7 


Calzetti law 


1067 


12 ±2 


8 


X-Ray luminosity 


24 


4 


9 


X-ray luminosity 




6.31 


10 


UV slope 


906 


5.7-18 


This Work 



References.— (1) Madau et al. 1996; (2)Meurer et al. 1997; (3) Pettini et al. 1998; (4) Meurer, 
Heckman, & Calzetti 1999; (5) Steidel et al. 1999; (6) Adelberger & Steidel 2000; (7) Teplitz et al. 2000; 
(8) Massarotti, Iovino, & Buzzoni 2001b; (9) Seibert, Heckman, & Meurer 2002; (10) Nandra et al. 2002 



by applying our lower attenuation limit (5.7) to 
the SFR and the upper one by applying our up- 
per attenuation limit of 18.6 to the observed SFR. 
If the upper limit applies, it implies that the vi- 
sual/UV surveys detect all the star-formation at 
this redshift, but, if the lower limit applies, it in- 
dicates a population of dust enshrouded galaxies 
which are not seen. 

Low-luminosity galaxies at high-z which may or 
may not be very dusty are probably missed by cur- 
rent UV/optical surveys. Hence they would also 
be missed in the correlation between the attenua- 
tion and luminosity of the systems at high-z. But 
if the dust formation in high-z systems is related to 
the evolution of massive stars and not of low-mass 
stars as it is currently in the Milky Way, then these 
galaxies will be intrinsically very luminous and be 
probably be detected by the UV/optical surveys. 

7. CONCLUSIONS 

We have presented our analysis of the dust in 
Lyman break galaxies at z <~ 3. Our sample con- 
sisted of 906 sources with spectroscopic redshifts, 
not identified as AGN or QSOs. Our main conclu- 
sions are as follows: 



1. The wide range of the G — 1Z colors of the 
LBGs, with most of them having similar 
apparent magnitudes suggests a major in- 
fluence of dust reddening. Possible varia- 
tion in chronological ages of the star- forming 
episode seems to be insufficient to account 
for this wide range of colors. 

2. We have used radiative transfer models in- 
corporating SMC-like dust in a shell geome- 
try (consistent with the Calzetti law) which 
also includes the dependence of the atten- 
uation function on the dust optical depths. 
The LBGs can be represented by an inner 
dust free region with high SFRs surrounded 
by a shell of clumpy dust. 

3. The observed G — 1Z colors appear to be 
incompatible with the presence of MW-like 
dust in these systems. The dust in LBGs 
is likely to be SMC-like, as found in other 
star-bursting systems. 

4. Complete resolution of the age-reddening de- 
generacy is not possible with a single ob- 
served rest-frame UV-color. However, the 
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Fig. 12. — Star formation rate history of the uni- 
verse. The lines are from Chary & Elbaz (2001), 
the solid line is for models that allow for lumi- 
nosity evolution, the dashed line (dotted in the 
electronic version) is the one for density evolution 
and the dotted line (dashed line in the electronic 
version) allows for both luminosity and the den- 
sity evolution. The crosses with the dotted error 
bars (colored blue in the electronic version) are the 
observed SFRs (Steidel et al. 1999), the solid er- 
ror bars (colored red in the electronic version) are 
our attenuation corrected values, with the lower 
limit obtained by applying the Case II correction 
of a factor of 5.7 and the upper limit obtained by 
applying the Case I correction factor of 18.4. 



dust attenuation can be constrained by con- 
sidering limiting cases where dust formation 
follows evolutionary timescales of massive 
stars (case I) and where the dust timescales 
follow evolution of intermediate mass (2.5 — 
3 Mq) stars (case II). The sources of dust 
production in case I are most likely to be 
SNe II. In case II the dust is most likely to 
be produced during the pre-planetary neb- 
ula stages of AGB stars. 

5. We performed an independent test of our 
reddening/attenuation models with limited 
set of LBGs for which V — H colors are avail- 
able. This test points toward an identical 
range of UV attenuations as found with the 
larger G — 1Z color set. 

6. Case I requires IMFs with upper mass limits 
increasing up to 200 M at 3.5 < z < 4.0 
and initial mctallicity Z = 0.01 Zq in order 
to match the bluest observed colors found 
among the LBGs. 

7. Both limiting cases I and II result in atten- 
uation factors which exhibit a pronounced 
dependence on the luminosities of the LBGs: 

a 160 o (Case I) oc (L 1600 A /L )°- 9O±a55 

a 1600 (Case II) cc (L 1600 A /L )°' 97±0 - 49 

8. We do not find any z-dependence in the dust 
properties of the LBGs up to z ~ 3.5, i.e. 
galaxies of similar luminosity are attenuated 
by similar amounts over the range 2 < z < 4. 
This suggests that the dust present in these 
systems is a direct result of the currently on- 
going star formation episode. 

9. The LBGs suffer a wide range of attenua- 
tions from zero to five mag; examined under 
both the limiting cases. 

10. Our luminosity weighted average attenua- 
tion factors for the rest-frame luminosity at 
1600 A for the LBG sample are 18.5 (case I) 
and 5.7 (case II). 

1 1 . When applied to the observed star formation 
rates at z <~ 3 and z <~ 4, the attenuation 
corrected star formation rates are consistent 
with model predictions based upon galaxy 
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counts and the cosmic IR background. This 
implies that the bulk of the star formation 
during the epoch 2 < z < 4 is directly 
detectable by current optical surveys, but 
requires substantial correction for internal 
dust attenuation. 
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